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The electronic states of the FeAs plane in iron-based superconductors are investigated on 
the basis of the two-dimensional 16-band d-p model, where the tight-binding parameters are 
determined so as to fit the band structure obtained by the density functional calculation for 
LaFeAsO. The model includes the Coulomb interaction on a Fe site: the intra- and inter-orbital 
direct terms U and U' , the exchange coupling J and the pair-transfer J'. Within the random 
phase approximation (RPA), we discuss the pairing symmetry of possible superconducting 
states including s-wave and d-wave pairing on the U'-J plane. 
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The newly discovered iron-based superconductors 1 ) 
RFeAsOi-zFz (R=Rare Earth) with a transition tem- 
perature up to T c — 55K 2 ) have attracted much at- 
tention. The mechanism of the superconductivity is one 
of the most significant issues. Several theoretical ap- 
proaches have been done on simplified multi-orbital Hub- 
bard models within weak coupling approaches. 3,4 ^ The 
details of the band structure and the Fermi surface are 
crucial for determining the pairing symmetry. Therefore, 
we employ a realistic model which includes both the Fe 
3d orbitals and the As 4p orbitals, so called d-p model. 

We perform the density functional calculation for 
LaFeAsO with the generalized gradient approximation 
of Perdew, Burke and Ernzerhof 5 -* by using the WIEN2k 
package, 6 ^ where the lattice parameters (a = 4.03268A, 
c = 8.74111A) and the internal coordinates (zj, a — 
0.14134, zas = 0.65166) are experimentally deter- 
mined. 7 -* Considering that there are two distinct Fe and 
As sites in the crystallographic unit cell, we then de- 
rive the two-dimensional 16-band d-p model, 8, 9 J where 



3d orbitals (d 3z 2_ r 2, d x 



d xy , d yz , d zx ) of two Fe 



atoms (Fei=A, Fe2=-B) and 4p orbitals (p x , p y , p z ) of 
two As atoms are explicitly included. We note that x, y 
axes are rotated by 45 degrees from the direction along 
Fe-Fe bonds. The model is given by the following Hamil- 
tonian, 

H = H Q + H inU (1) 
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where duo- is the annihilation operator for Fe-3d electrons 
with spin a in the orbital £ at the site i and Pi ma is 
the annihilation operator for As-4p electrons with spin 
a in the orbital m at the site i. In eq. (2), the transfer 
integrals tff j e e ,, t pp j m m ,, tf p e m and the atomic energies 
ef, e^i are determined so as to fit both the energy and the 
weights of orbitals for each band obtained from the tight- 
binding approximation (d-p model) to those from the 
density functional calculation. The doping x corresponds 
to the number of electrons per unit cell n = 24 + 2x in 
the present model. 

In eq. (1), Hi nt is the on-site Coulomb interaction 
in Fe-3d orbitals and includes the intra-orbital (inter- 
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Fig. 1. Fermi surface obtained from the d-p model eq. (2) for 
x = 0.1 
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Fig. 2. Several components of the spin susceptibility obtained 



a,/3 



from the RPA X e iz 2 , , 
J = J' = 0.25eV, T = 0.02eV and x 



(q) for U = 1.5eV, V 



l.OeV, 



0.1. 



Fig. 3. The gap function A^ 



(k) for extended 



-y* — y* 

(left) with A = 0.976 and for d xy -wa,ve (right) with A = 0.836 
for the same parameters in Fig. 2. The solid line represents the 
node of the gap function. 



orbital) Coulomb interaction U (V), the Hund's rule 
coupling J and the pair-transfer interaction J'. In the 
weak coupling regime, the superconducting gap equation 
is given by 

*a#ck)=4e E EE 

k' e 1 e 2 £ 3 i 4 a', 13' li,v 

x u^' )/1 (-k')«? liP (-k / )*< ) „(k , )«f !liI/ (k')* (3) 

where /i, v (=1-16) are band indexes, a, (3 (=A, B) repre- 
sent two Fe sites, u"„(k) is the eigenvector which diago- 
nalizes Hq, £k.p is the energy of band /i with wave vector 
k, /(e) is the Fermi distribution function. Here, we ap- 
proximate the effective pairing interaction Vj>™^£ ^ (c|) 
within the RPA 10 ' and numerically solve the equation 
(3) to obtain the gap function A"^(k) with eigenvalue A. 
At T = T c , the largest eigenvalue A becomes unity. We 
use 32 x 32 k points for numerical calculations. 

Figure 1 shows the Fermi surface for the d-p model 
at x = 0.1, where we can see nearly circular hole pockets 
around the T point (k = (0, 0)) and elliptical electron 
pockets around the M point (k = (it, tt)). The spin sus- 
ceptibility X^ 1 " 2 '^ 3 ^ 4 ( t l) is given by 50 x 50 matrix and is 
calculated within the RPA as shown in Fig. 2. The peak 
around the T point is due to the effect of the nesting be- 
tween the electron (hole) pockets, while the peak around 
the M point is due to that between the hole pockets and 
the electron pockets. We note that the x 2 —y 2 component 
becomes dominant around the M point (see Fig. 2). 

As the enhanced spin susceptibility contributes to 
the effective interaction within the RPA, 10 ) the gap func- 
tion for the x 2 —y 2 component A^_ y2 aJ 2_ y 2(k) becomes 
dominant. The obtained gap functions with two largest 
eigenvalues A are shown in Fig. 3. For U = 1.5eV, 
U' = l.OeV, J = J' = 0.25cV, T = 0.02eV and x = 0.1, 
the gap function with the largest A is extended s-wave 



symmetry and changes the sign between the hole pockets 
and the electron pockets, while that with second largest A 
is (i^y-wave symmetry, as shown in Fig. 3. On the other 
hand, for U = 1.68eV, U' = 1.4eV, J — J' — O.UeV, 
T = 0.02eV and x = 0.1, the gap function with largest A 
is d^y-wave (A = 0.978), while that with second largest 
A is extended s-wave (A = 0.967). The detailed phase 
diagram of the pairing symmetry will be shown in the 
subsequent paper. 9 ) 

In summary, we investigated the pairing symmetry 
of the two-dimensional 16-band d-p model by using the 
RPA. For a larger value of J/U', the most favorable pair- 
ing is extended s-wave symmetry whose order parame- 
ter changes its sign between the hole pockets and the 
electron pockets, while for a smaller value of J/U', it 
is c^y-wave symmetry. According to the recent experi- 
ment of very weak T c -suppression by Co-impurities, n ) 
we suppose that the c^y-wave pairing is suppressed by 
pair breaking effect and the extended s-wave paring is 
realized in real materials. 
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